Quantum quenches in a spinor condensate 
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We discuss the ordering of a spin-1 condensate when quenched from its paramagnetic phase to 
its ferromagnetic phase by reducing magnetic field. We first elucidate the nature of the equilibrium 
quantum phase transition. Quenching rapidly through this transition reveals XY ordering either at 
a specific wavevector, or the 'light-cone' correlations familiar from relativistic theories, depending on 
the endpoint of the quench. For a quench proceeding at a finite rate the ordering scale is governed 
by the Kibble-Zurek mechanism. The creation of vortices through growth of the magnetization 
fluctuations is also discussed. The long time dynamics again depends on the endpoint, conserving 
the order parameter in zero field, but not at finite field, with differing exponents for the coarsening 
of magnetic order. The results are discussed in the light of a recent experiment by Sadler et al. 
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How does a many-particle system undergo condensa- 
tion into an ordered state? This question is central to 
a number of disparate areas of physics, from condensed 
matter to cosmology P, 0] . Often we are interested in the 
processes determining the formation of ordered domains 
and topological defects. The usual approach is to study 
the coupled dynamics of the collective (or hydrodynamic) 
degrees of freedom, including the order parameter and 
any conserved quantities. Thus the dynamics is highly 
constrained by the presence or absence of conservation 
laws, with dramatic differences in the resulting time evo- 
lution of correlations. Such considerations will be im- 
portant when we quench into an ordered state at zero 
temperature through a quantum phase transition [3, . 

Such a possibility was explored in a recent experiment 
that studied ferromagnetic ordering in a Bose-Einstein 
condensate of ^^Rb atoms following a sudden reduction 
in magnetic field ^ . Cold atomic gases represent an ex- 
citing new prospect for the investigation of such quantum 
quenches. As we will show, they represent a far closer 
analog of relativistic theories than the condensed matter 
systems suggested for 'laboratory cosmology' by Zurek 
and reviewed in Ref. Q 

Earlier, mostly numerical work 

0, BSE El has fo- 
cused on treating the creation of spin domains in conden- 
sates as a property of a classical dynamical system. Our 
goal in this Letter, on the other hand, is to first explain 
the character of the equilibrium quantum phase transi- 
tion, and then to discuss the associated dynamics as a 
problem of phase ordering. In particular, this will lead 
us to carefully distinguish different quenches in terms of 
the conservation laws obeyed, and the resulting dynam- 
ics of topological defects (vortices in the magnetization) . 
These are features not present in the model systems dis- 
cussed in Refs. 0, [^, for instance, and their treatment 
requires the introduction of some novel theoretical ideas 
that should be of broad applicability. 

The existence of an ordering transition in the ferro- 




FIG. 1: Zero temperature phase diagram of a spin-1 conden- 
sate in terms of the linear and quadratic Zeeman energies. 
The shaded area corresponds to the region of XY ordering. 
The experiment of Ref. [3| involved a quench through the spe- 
cial point at p = (red dot) from the paramagnetic region. 



magnetic spin-1 Bose gas is readily understood on the 
basis of a variational Gross-Pitaevskii calculation Il2|. 



The second-quantized Hamiltonian is [13|, [14 1 



H = / drHo+Hir, 
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Hint = -j^4'ln4'ln' [cO^ranSm'n' + C2F nin ' 'P m' n'] (f'n''^n 

The index m = —1,0, 1 gives the z-component of total 
spin of the corresponding state, and we have set h and 
the atomic mass to unity. are the spin-1 matri- 

ces for i = x,y, z. E^^ is the Zeeman energy of the m- 
component, defined below. For a ferromagnetic system 
the spin interaction parameter C2 is negative. We imple- 
ment the Gross-Pitaevskii approximation by treating the 
0„ as c-numbers tpm and writing ifm — V^Xm in terms 
of a normalized spinor Xm- The energy per particle is 
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Con C2n 



(IXiP - IX-il')' + 2 (IxoXiP + IXoX-iP + Xo'xiX-i + XoXtxli) + 9 [IXiP + IX-iP] - P [IXiP - Ix-i 



(2) 



The last two terms originate from the hnear and 
quadratic Zeeman energies, the general case for spin-1: 
E"^ = —pm + qm? . In fact the coefficient p in the energy 
is the sum of the linear Zeeman term and a Lagrange mul- 
tiplier enforcing conservation of . We should minimize 
Eq. ^ and then use the true value of the magnetization 
{Fz) to fix p. Thus with zero total magnetization we will 
have p = 0. In this case it is straightforward to see that 
for q > qo = 2\c2\n the spinor state x' = (0, 1,0) mini- 
mizes the energy, while for q < qo the ±1 states become 
populated, leading to a transverse magnetization density 
/x(r) ^ r(r) +i/nr), where ^(r) ^ 0L(r)i^^>„(r). 
In the general case p ^ [l5| 



showing that {f±) is nonzero between the lines \p\ — q 
and \p\ = \/qiq — qo) (see Fig. [T]). Thus the mean-field 
calculation predicts a symmetry-breaking quantum phase 
transition. For p ^ 0, there is also a perpendicular mag- 
netization in this region 



p (p2 + q^q _ q^) 



qoq 



and the ordered phase is a canted XY ferromagnet [iq, 
while the /_l = phase is an XY paramagnet (finite 
transverse susceptibility). 

We now ask in more detail what kind of quantum phase 
transition we are dealing with. This is more than a for- 
mal question, as the dynamics of the order parameter 
at the transition will be crucial in determining the be- 
havior following a quench. We start by considering the 
Bogoliubov theory of the paramagnetic phase. Shifting 
the fields 0^ in Eq. ^ hy ip^ = (0, V", 0) we find that 
in the quadratic part of the Hamiltonian the states 
decouple from the 0q state to give 



Hb^ 4>l^^ [ek-pm + q + C2n] (j)^^^ 

A:,m— ±1 

k 

{sk = fc^/2) This is readily diagonalized by Bogoliubov 
transformation to yield (except for a constant) 



(5) 



where Es^±{k) = Es{k)^p, with the spin wave dispersion 
-^s (^) ^ (sfe + q) i^k + q + 2c2n). One of these disper- 
sions passes through zero when \p\ — pc = \/ q{q — qo), 
the same instability of the paramagnetic phase that we 
found before. Except at p = the transition to the or- 
dered state proceeds by filling of either the 'particle' or 
'hole' band in Eq. ([5]), and involves a change in Fz- Since 
this is conserved, it is impossible to cross this transition 
without contact to a reservoir of magnetization (unless 
T 7^ 0). When p = 0, on the other hand, the transition 
occurs through closing of the bandgap and the longitu- 
dinal magnetization remains zero. 

To make a connection to the general theory of quantum 
phase transitions, we rewrite the Hamiltonian Eq. ^ 
using the (complex) canonical coordinates and conjugate 
momenta 
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+ blk) . ^k = ^VeM (4 - b-k) 



In this way we get (dropping the momentum sum) 



-ffs = ^ (tt^ -I- ipz) (tt 



ipz^ 



(6) 



Eq. ^ is recognized as the Hamiltonian of a two- 
dimensional particle in a uniform perpendicular magnetic 
field and harmonic oscillator potential. 

How do we interpret the field z? The Fourier modes 
of the transverse magnetization density /±(r) may be 
written in terms of as 



= \/2n{ek + q)zl -I 



4+k,0^l-l 



where the dotted lines denote terms higher order in the 
quasiparticle operators. At low k the two are simply 
proportional, as one would hope. 

Below the transition, the higher order terms dropped 
from the Hamiltonian Eq. ([4]) are required to saturate 
the growth of /_L(r). Close enough to the transition a 
quartic term is sufficient. Obtaining this term within 
the Bogoliubov theory is slightly subtle as it involves 
the partial cancellation of the 'direct' quartic interac- 
tion neglected in Eq. ([4]) against the interaction induced 
by phonons [17I] . The most relevant term may be found 
more simply, aside from a small renormalization, from 
the observation that it is responsible for the square root 
singularity in Eq. ([3]) in the mean-field approximation. 
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We write the final result as an effective action valid near 
the transition, and for length scales ^ |27[. 



1 



S^ff = I drdt- [z^z - cl\7z''Vz - 



[Pc 



-ipz^ z - 



p2) ^t^l 



^1 ,4 



(7) 



C2 = (1 — ■^)c2. To obtain the quadratic part of Eq. (O 
from Eq. (0) we have approximated the spectrum as 



Pc 



= q-qo/2 



(8) 



The precise form of the anharmonic terms in Eq. ([7]) will 
not be important in the following development. 

^eff is identical to the effective theory describing 
the superfluid-insulator transition in the Bose-Hubbard 
model. As in that problem, the point p = is identified 
as a special point where the transition, instead of being of 
the bose condensation type, lies in the universality class 
of the {d + l)-dimensional XY model It is relevant 
to ask whether the deviations from mean field critical be- 
havior implied by this identification will be seen in exper- 
iment. For a two-dimensional condensate a standard cal- 
culation gives the Ginzburg criterion for the breakdown 

1/2 

of mean field behavior O 4^ f > 1 where 
Lj_ is the transverse dimension. Since for the system in 
Ref. 5] the prefactor is of order lO^"', the mean field 
theory is an excellent approximation. 

With this mind, we now proceed to describe the evo- 
lution of a system that is quenched suddenly through 
the transition at p = from an initial value qi > qo for 
t < to (7/ < go at t > 0, as in Ref. There is a 
band of unstable modes with E'^{k,qf) < 0, the occu- 
pancy of which begins to grow exponentially, as they are 
populated with pairs of atoms scattering from the m = 
state. The quadratic Hamiltonian Eq. (O describes this 
process adequately until the populations are such that 
the anharmonic interactions between the modes - such 
as the last term of Eq. [7]- become important. Writing 
the Hamiltonian for t > in terms of the fields Zf, defined 
at g = (7i, and with w| = —E^(k, qf), we find the solution 
of the Heisenberg equations of motion 



= Zj,(0) coshwfei + 7r^(0)wj, sinhwfci 



The calculation of the correlation function of the trans- 
verse magnetization is then staightforward 

(/±feW/LfeW> = 2n(efe +qO \cosh^u;kt{zl{0)z^{0)) 



£k + qi 



sinh^Wfet(^^(o)^jO)) 



revealing the exponential growth of the magnetization. 
The initial fluctuations of the oscillator modes are [2^ 

(zt(0)z,(0)) = {Es(k,q,))-\ (7rt(0)7r,(0)) = E,{k,q,). 



These general formulae are valid for any instantaneous 
quench (29| . In the following we will make the simplifi- 
cation of taking qi ^ go (shot noise limit) , which gives 



cosh LUkt - 



£fc + qf 



sinh LUkt 



Now we wish to focus on two particular values of qf to il- 
lustrate the different possible classes of behavior. If qf = 
the spectrum of unstable modes is uj'^ — Sk (go ~ ^k), 
which has a maximum ai k — The correlation func- 
tion is therefore dominated by the fluctuations on this 
scale that grow at a rate go. Taking into account only 
the unstable modes, we find for the asymptotic behavior 
of the real space correlations 




MVqol^ - r'l)e 



qot 



(/^(r,t)/|(r',t)) 



for qot ^ 1, r/t <C Cg. The Bessel function is an angular 

1 /2 

average of plane waves of wavevector q^ . The result is 
a growing random spin texture of typical scale go ^^"^j as 
observed in Ref. (Hj . Note that the vanishing of the mode 
growth rate at zero wavevector is a consequence of the 
conservation of all three spin components in zero field. 

Very different behavior results if g/ is only just below 
go, so that Pc/c1 << 1. In this case the spectrum of 
unstable modes reflects the relativistic form of Eq. ([8|) 

Lui - ci {ki - e) , 

with k1 = —p^/Cg the 'Compton wavevector'. In this 
case we get the asymptotic behavior 



(/^(r,t)/|(r',t)) 



ngo 



1 



ATrL± Cgkct 



(9) 

valid when the exponent is large. The correlation func- 
tion Eq. Q displays a striking growth of correlations 
along a 'light cone' originating at a point halfway be- 
tween r and r' and propagating at the spin wave velocity 
Cs- This is a familiar feature of spinodal instabilities in 
relativistic theories [l^. The crossover between these 
two types of behavior occurs at the value g/ = go/2, 
where the maximum in the spectrum of unstable modes 
fcmax = V90 - 2g goes to zero. 

Next we discuss what happens if the quench is not in- 
stantaneous, but rather crosses the transition in some fi- 
nite time. For concreteness we take q(t) = go (1 — t/Tg), 
where tq measures the duration of the quench, and the 
transition is crossed at < = 0. The result may be obtained 
exactly in terms of Airy functions, but the following in- 
tegral representation is more useful 



(/^(r,i)/l(r',t)) 



nqo 



2citKzL± 



T{t/tKz,\r~r'\/{cstKz)) 



r3/2 



exp 



Au 12 
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where we have introduced t^z = {tq/QoY^^ ■ This ex- 
pression can then be evaluated in the saddle-point ap- 
proximation. At |r — r'l = we get an exponential factor 

exp ^1 (i/iKz)'^^^y At this point we have to invoke for 
the first time the effect of the anharmonic interactions 
between modes. Very crudely, their effect is to cut-off 
the exponential growth of the magnetization. We shall 
not try to discuss this process in detail, but the key point 
is that it occurs at a time cx Ikz, where the constant of 
proportionality may contain t^z, but only logarithmi- 
cally. Thus we readily see that the associated scale is 

1 /3 

£,{tKz) oc CgtKZ = Cs {tq/Qx)) ■ This result is consistent 
with the general arguments of Kibble and Zurck, imply- 
ing a domain size scaling as (tq) , with mean field 
values z = 1 and v = \j1 for the dynamic and correla- 
tion length exponents [6[ . 

The growth of the transverse magnetization is asso- 
ciated with the appearance of vortices. As the popu- 
lation of the unstable modes becomes large, the field 
/j^(r) can be treated as an effectively classical Gaus- 
sian stochastic variable, with variance given by the cor- 
relation functions calculated above [20]. Then the den- 
sity of vortices can be estimated using the Halperin-Liu- 
Mazenko formula to calculate the density of zeroes of 
this classical field nv{t) = -^dlg" {r,t)\r^Q [HI, E^, 
where 5(|r-r;|,i) ^ (/^(r, i)/l(r', i))/(/^(0, t)/|(0, t)) 
is the normalized correlation function. For quenches 
to qf < |c2|n it is immediately clear that the den- 
sity is determined by fc„iax, as the spectrum of fluctu- 
ations is essentially monochromatic at late times, and 
k'^^^/An. In this case the vortices have a 
core size of the same order as the scale of the mag- 
netic order. For a quench to just below the transi- 
tion, on the other hand, the asymptote Eq. ([9|) gives 
nv{t) ^ j^-^- This behavior continues until the growth 
is saturated by the anharmonic terms, which happens 



when l/i 



2n- (<Zo — q) /c2- Finally, in the case of the 



finite time quench, we have ny cc C^^(iKz)- 

In closing, we discuss the long-time behavior of the 
system, once the transverse magnetization is comparable 
to its equilibrium value. This regime is characterized by 
the growth of the characteristic ordering scale and the 
annihilation of topological defects, usually called coars- 
ening. We distinguish two universality classes depending 
on whether or not the order parameter is conserved [ij . In 
the first case the domain size increases as ^(i) oc t^/^ [30|, 
while in the second a t^^^ law is obeyed. In our system 
these two cases correspond to a final value ^ or 
Qf = respectively. Note that for the shallow quench we 
found the t^^^ behavior already at the linear level. 

A further complication is that coarsening is usually 
studied using models of dissipative dynamics, where en- 
ergy is not conserved. On the other hand, coarse-graining 
of a purely Hamiltonian system can give rise to such dy- 
namics, at the expense of introducing a conserved energy 



density to which the order parameter is coupled [23[ . For 
the case of a real scalar non-conserved order parameter, 
Hamiltonian coarsening was examined Ref. [2J] , with the 
conclusion that the t^^^ law was preserved (this case cor- 
responds to Model C in the classification of Ref. |i23i]). 
On the other hand, Ref. [2^ studied coarsening in the 
Gross-Pitaevskii equation, where energy and addition- 
ally particle number are conserved (Model F), and found 
results consistent with ^(t) cx t. The dynamics described 
by Eq. ^ corresponds to Model F, except at the particle- 
hole symmetric point p = that has been our main con- 
cern, a special case called Model E. 

I would like to thank John Chalkcr, Joel Moore, Sab- 
rina Leslie, and Julien Kockelkoren for useful discussions. 
After this work was finished, the preprint Ref. [26\ ap- 
peared, where the same experiment is discussed. 
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factors coth {Es{k) /2kBT). [30] For a complex order parameter in d = 2 the presence of 

[29] Note that the hnear Zeeman term leads to a trivial Lar- vortices implies ^(t) oc (t/lnt)^''^. 

mor precession of the magnetization that is not included. 



